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Abstract 

Anderson & Schubert (2007, Science, 317, 1384) proposed that Saturn's rotation period can be as- 
certained by minimizing the dynamic heights of the 100 mbar isosurface with respect to the geoid; they 
derived a rotation period of lOh 32m 35s. We investigate the same approach for Jupiter to see if the 
Jovian rotation period is predicted by minimizing the dynamical heights of its isobaric (1 bar pressure 
Oh level) surface using zonal wind data. A rotation period of 9h 54m 29s is found. Further, we investigate 

the minimization method by fitting Pioneer and Voyager occultation radii for both Jupiter and Saturn. 
Rotation periods of 9h 55m 30s and lOh 32m 35s are found to minimize the dynamical heights for Jupiter 
Oh an d Saturn, respectively. Though there is no dynamical principle requiring the minimization of the dy- 

namical heights of an isobaric surface, the successful application of the method to Jupiter lends support 

5-h to its relevance for Saturn. 

-I— > 

We derive Jupiter and Saturn rotation periods using equilibrium theory in which the solid-body rotation 
period (no winds) that gives the observed equatorial and polar radii at the 100 mbar level is found. 
Rotation periods of 9h 55m 20s and lOh 31m 49s are found for Jupiter and Saturn, respectively. We 
show that both Jupiter's and Saturn's shapes can be derived using solid-body rotation, suggesting that 
zonal winds have a minor effect on the planetary shape for both planets. 

The agreement in the values of Saturn's rotation period predicted by the different approaches supports 
^j- the conclusion that the planet's period of rotation is about lOh 32m. 
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1 Introduction 



Knowledge of a planet's rotation rate is fundamental to understanding its internal structure and atmospheric 
dynamics. The atmospheres of Jupiter and Saturn have strong zonal winds with equatorial speeds reaching ~ 
100 m s _1 on Jupiter and ~ 400 m s _1 on Saturn. However, these zonal wind velocities are based on assumed 
values of the rotation periods of the planets, 9h 55m 29s for Jupiter, and lOh 39m 22s for Saturn (Lindal, 
1992). Jupiter's rotation period is assumed to be identical to the rotation period of its tilted magnetic field, 
which has not changed in many decades (Porco et al., 2003; Smith & Hunt, 1976, Riddle & Warwick, 1976; 
Higgins ct al., 1996). The lOh 39m 22s rotation period of Saturn, on which the 400 m s _1 equatorial zonal 
wind velocity is based, derives from the Voyager era periodicity in Saturn's kilometric radiation (Ingersoll 
& Pollard, 1982; Dessler, 1983). That period is now known to be variable and Saturn's rotation rate is 
therefore uncertain (Gurnett et al., 2007). Accordingly, the atmospheric zonal wind velocities with respect 
to the underlying rotating planet are not known for Saturn. 

Anderson & Schubert (2007, hereafter paper I) used the improved gravitational data measured by Cassini 
(Jacobson et al., 2006) to infer the properties of Saturn's interior. Since the rotation rate of Saturn is 
unknown internal models were presented for rotation periods between lOh 32m 35s and lOh 41m 35s. The 
interior structure was modeled by applying the 'theory of figures' (Zharkov & Trubitsyn, 1978) using a 
6th order (in normalized radius) polynomial density function, assuming a solid-body rotation. Anderson & 
Schubert (2007) suggested that Saturn's rotation period could be constrained by minimizing the dynamical 
height variations with respect to the geoid of the 100 mbar isosurface of Saturn's atmosphere. The shape 
of Saturn's 100 mbar isosurface is known, but the shape of the geoid is not because it depends on Saturn's 
unknown rotation rate. A minimization of the dynamic heights of the 100 mbar isosurface was found for a 
rotation period of lOh 32m 35s. This rotation period enables computation of an interior model that best 
fits the measured gravitational moments, and is consistent with the lower bound rotation period from all 
available magnetic and SKR data (Gurnett et al., 2007). Figure 1 shows the latitude-dependent dynamic 
heights of the 100 mbar isobaric surface, measured relative to a reference geoid with rotation periods of lOh 
32m 35s, lOh 35m 35s, and lOh 38m 35s (paper I). 

Smith ct al. (1982) commented that a more rapid Saturnian rotation rate (compared with the Voyager 
rate) would result in lower wind velocities with respect to the solid-body rotation, and that the dynamical 
heights are minimized for a rotation period between lOh 31m and lOh 32m. A shorter rotation period 
for Saturn's deep interior has also been suggested by Allison & Stone (1983) who noted that a rotation 
period shorter than the Voyager radio rotation period by about 1% could better explain the energetics 
of the Saturnian atmosphere. A shorter rotation period for the interior produces a Jovian-like Saturnian 
wind system with sensible equatorial wind speeds and both eastward and westward jets at higher latitudes. 
However, at that time, and during the last few decades there was no reason to question Voyager's rotation 
period, and the idea of a shorter rotation period for Saturn was not pursued. 

Cassini has recently measured the SKR period to be lOh 47m 6s (Gurnett et al. , 2007) , about 8 minutes 
longer than the Voyager period. In addition, during Cassini's orbit around Saturn this period was found to 
be changing with time (Gurnett et al. , 2007) . It is therefore clear that these measurements do not represent 
the rotation period of Saturn's deep interior. The good fit of the gravitational data for an interior model 
with a faster rotation rate than the Voyager rate, as shown in paper I, and the realization that Saturn's 
rotation period is not given by the SKR period, call for new efforts to confirm the proposed faster rotation 
of Saturn. 

In this paper we present different methods for inferring the rotation periods of the giant planets. Ap- 
plying the methods to Jupiter and taking Jupiter's rotation period as unknown, we reconstruct Jupiter's 
SKR rotation period with an error of up to 1 minute. We find that the dynamical heights minimization 
method (paper I) works successfully for Jupiter. However, there is no dynamical principle requiring the 
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minimization of the dynamical heights of the isobaric surface and the method does not work for terrestrial 
planet atmospheres. The minimization method docs work for Jupiter whose atmosphere has both eastward 
and westward jets. The successful application of the method to Jupiter lends support to its relevance to 
Saturn, for which a rotation period of ^10h 32m is derived. 



2 Minimization of the Dynamical Heights of a Pressure Isosurface 
Using Wind Data 



Jupiter's physical shape, i.e., the radius of its isobaric surface (here taken as the 1 bar isosurface) as a 
function of latitude, can be represented by a combination of its reference geoid and the dynamical heights 
of the isobaric surface relative to the geoid. The reference geoid is the spheroid whose shape is determined 
by equilibrium theory for a given rotation period (Zharkov & Trubitsyn, 1978). Using zonal wind data 
with respect to an assumed uniform rotation period, the radius of the corresponding potential surface as a 
function of latitude can be computed (Lindal et al., 1985; French et al., 1998). The isobaric surface does not 
coincide with the geoid because of the dynamical effects of zonal winds. The dynamical heights with respect 
to the reference geoid are given by (Lindal et al., 1985) 



where g is the average of the magnitude of the acceleration of the gravity (the vector acceleration of gravity 
varies both radially and latitudinally along a local vertical between the isobaric surface and the reference 
geoid), V is a small angle that gives the difference between the planetocentric latitude cf) and the planeto- 
graphic latitude <fi', ip = <p' — <p, r is the distance to the center of the planet, and Vw is the zonal wind 
velocity. Parameters with the subscript 'ref ' represent quantities calculated for the reference geoid (sec Lin- 
dal et al., 1985 for further details). For simplicity, and because vertical profiles of zonal wind speeds are 
unknown, the measured zonal wind velocities are taken to be independent of altitude (radial component). 

For Vw we use zonal wind data obtained from Hubble Space Telescope images during 1995 to 2000 
(Garcia-Melendo & Sanchez-Lavega, 2001) at three different wavelengths (410, 892, 953 nm) representative 
of Jupiter's cloud top levels. Figure 2 shows the wind velocities for Jupiter as computed by Garcia-Melendo 
& Sanchez-Lavega (2001). These wind data are found to be in good agreement with previous zonal wind 
measurements for Jupiter (Limaye, 1986). 

We calculate Jupiter's reference geoid for its measured rotation period of 9h 55m 29.7s (Lindal, 1992). 
The reference geoid is defined as the level surface of equal effective potential 



The level surface of constant effective potential U is computed to the fifth order in the smallncss parameter 

2 3 

q = ^ff, where a is the equatorial radius of the geoid, GM is its mass multiplied by the gravitational 
constant, 9 is the co- latitude, and u is the angular velocity given by 2n/P, with P the rotation period. 
The polar radius is fixed at the value 66,854 km, the polar radius of the 1 bar isosurface. The value of the 
polar radius is taken as the average of the two polar radii at that pressure level. The northern and southern 
polar radii differ from each other due to the difference in the zonal wind strength at the two hemispheres 
(e.g., Ingersoll, 1970). The harmonic coefficients Jm are obtained from the measured values for a reference 
equatorial radius of 71,492 km (Jacobson, 2003). The values of Ji n from Jacobson (2003) in units of 10~ 6 




(1) 





are J 2 = 14696.43 ± 0.21, J 4 



587.14 ± 1.68, J 6 = 34.25 ± 5.22. 
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The measured gravitational coefficients J 2n at a reference radius of 71,492 km arc of order q n , and can 
be replaced in the effective potential function by a 2n j2n=(71,492 km) 2 ™ J 2n (measured). Since the equatorial 
radius a of the reference geoid depends on the rotation period, so do the values of j 2n . 

The radius rasa function of cosO can be expanded by 

r{cos6) = a ^1 + ^ r n {cos6)q^j . (3) 

To find the reference geoid we substitute the expression for r (equation (3)) into equation (2) and expand 
the result in the smallness parameter q. Then a zero function can be defined by, 

F(cos6) = U{cos6) - 17(0) (4) 

where, to the fifth order in q, 

u(0) = i + \ q + \m \m 2 + ^m 3 ^jW + §- 6 m q 5 . (5) 

The coefficients r n are then found recursively by setting the coefficients of q in the expanded equation (4) 
to zero, starting with order one. The results for the first three r n functions are, 

n = - l -{l + Zj2)cos6 2 

r 2 = - l - (2 + 3j 2 - 9j| - 15j 4 ) cosO 2 + 1 (6 + 6 (1 - 6j 2 ) 32 - 35j 4 ) cos6 4 
r 3 = - ^ (8 - 27 h (2 (1 - j 2 )j 2 - 5j 4 ) - 45j 4 + 105j 6 ) cosO 2 

+ 16 (6 ~ &ji + 3j2 (2 + 6j2 (3j2 ~ :) + 43j4) + 63je) C ° s9i 

- ^ (24 + 35.74 + 18.72 (2 + J 2 + 21j| + 35j 4 ) + 231 j 6 ) co S 6 (6) 

The expressions are considerably more complicated for orders 4 and 5, but they can be evaluated similarly. 
The point in going to fifth order is not to introduce measured values of Jg and Jio, they have not been detected 
so far; fifth order is used for purposes of representing the measured radii to a precision commensurate with 
their accuracy at the ~ 2 km level. 

This procedure provides an equation for a/r. We then determine a/b where b is the polar radius, by 
a/b = a/r(cos9 = 0). Since the polar radius is assumed to be known (66,854 km), we finally find the 
equatorial radius a, and the values of j2„ for an assumed rotation period ui (see paper I for further details). 

Once the reference geoid is calculated, Jupiter's physical shape can be derived. The physical shape 
represents the height of the 1 bar level with respect to the planet's center, and is constructed by adding 
Jupiter's dynamical heights h((f>) to the reference geoid corresponding to the measured rotation period of 9h 
55m 29.7s. We can now search for the rotation period which minimizes the dynamical heights of the isobaric 
surface (equation (1)). To do so, we construct reference geoids for four rotation periods: 9h 57m 29.7s, 9h 
55m 29.7s, 9h 54m 29.7s, and 9h 53m 29.7s. The altitudes of the 1 bar isosurface above the reference geoids 
(dynamical heights h(4>)) for different assumed rotation periods are found by subtracting the calculated 
geoids (with varying rotation periods) from the physical shape of Jupiter (determined from the dynamic 
heights with respect to the 9h 55m 29.7s geoid). 

Table 1 provides the physical parameters of the geoid for the four assumed rotation periods. The ampli- 
tudes of the 1 bar isosurface above the reference geoid for the four different rotation periods are shown in 
Figure 3. The dynamical heights for Jupiter's measured rotation period are found to be in agreement with 
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the ones calculated by Ingcrsoll (1970). The rotation periods investigated are 1 minute shorter and 2 minutes 
shorter and longer than the actual measured rotation period of Jupiter, because the inherent standard error 
in rotation period using the dynamical height minimization method was found to be about ±1 minute, for 
both Jupiter and Saturn. The standard error for the wind radii is a few kilometers resulting in an error 
in rotation period of up to 1 minute (paper I). As can be seen from the figure the dynamical heights are 
minimized for a rotation period of 9h 54m 29.7s, a rotation period which is only one minute faster than 
Jupiter's measured rotation period. For Saturn, the method results in a rotation period of lOh 32m 35s 
(paper I). 



2.1 Minimizing the Dynamical Height Difference Between the Equator and the 
Pole 

A simplified version of the dynamical height minimization method presented above minimizes the dynamical 
height difference between the planet's pole and equator. We take vy to be the measured cloud-tracked wind 
velocity, referenced to the Voyager radio period (Garda-Melendo & Sanchez-Lavega, 2001). The zonal wind 
measured with respect to a 'revised reference geoid' v r varies with the revised planetary rotation period as 
(Allison & Stone, 1983) 



v r = vy + 2ira ( I coscf) (7) 

\Ty T r J 

where a is the planet's equatorial radius, ry is the Voyager rotation period, and r r is the revised planetary 
rotation period. 

The dynamical heights with respect to the reference geoid (equation (1)) for small wind velocities and 
negligible small non-geostrophic contributions to the integral are given by (Smith ct al., 1982; Lindal et al., 
1985), 

h(4>) = j y vS i n( j)d(p } (8) 

9 J 4> 

where ujy is the Voyager angular velocity and g is the average acceleration of gravity. Then, by substituting 
equation (7) for the zonal wind into equation (8), the implied change in the dynamical height difference 
between the equator and the pole can be written as 

a 2 / 1 1 \ 

h r = hy + 2-KUJy— , (9) 



9 \ T V T r 

where hy is the Voyager dynamical height. Assuming that the dynamical heights are minimized (h r goes 
towards zero), equation (9) can be solved for St = r r — ry with the result 

x Shy , . 

5T = - TV Waf- (10) 

Substituting the values for Jupiter {hy <~ 12 km, g <~ 24 m s~ 2 , (ujya) 2 ~ 1.58xl0 8 m 2 s~ 2 , and ry — 
9.925h = 9h 55m 30s), we get St <~ -1 minute, essentially the same as the 1 minute shorter period derived 
earlier. 

Applying this approach to Saturn, with hy = 120 km, g ~ 10 m s~ 2 , (utya) 2 = 9.765 x 10 7 m 2 s~ 2 , and 
Ty = 10.6567h = lOh 39m 24s (Lindal et al., 1985), we get St <~ -7 minutes, a value in excellent agreement 
with the revised period of lOh 32m 35s derived in paper I. Although the uncertainty of these estimates is 
a few tens of seconds, we show that the simplified version of the dynamical height minimization method 
produces comparable results to the ones obtain above. 
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3 Fit to Pioneer and Voyager Occultation Radii 



The physical shape of Jupiter and Saturn can be derived independently of zonal wind data by using occul- 
tation radii. While in the previous section the physical shape was derived from the geostrophic equation, 
occultation radii provide the planetary shape independently of wind data. 

Six radio occultation radii of Jupiter's 100 mbar isosurface from the Pioneer and Voyager missions are 
given in Fig. 7 of Lindal et al. (1981). Using fifth-order expressions for the radius (equation (3)), we 
recreate the geoid of Lindal et al. (1981) by using the same gravitational field (Null, 1976) as used in their 
paper. The equatorial radius of this geoid is 71,541 km at the 100 mbar level, with a polar radius of 66,896 
km. To construct Jupiter's physical shape (occultation radii at their corresponding geocentric latitudes) 
independently of the underlying geoid, we add the heights of the 100 mbar isosurface with respect to the 
geoid by digitizing the right hand panel of Fig. 7 in Lindal et al. (1981) to two places past the decimal, well 
inside their error bars of ~ 2 km. 

We next search for a rotation period that can best fit the occultation radii (physical shape), i.e., a period 
that minimizes the heights between the calculated geoid (for varying rotation periods) and the occultation 
radii. We introduce three free parameters into the fitting function. The first is the equatorial radius a at the 
100 mbar level, the second is the period of rotation P, and the third is an offset AZ of the center of mass 
along the polar axis, as suggested by Lindal et al. (1981). This offset produces a better fit to the occultation 
radii. Both a and P are used to update the smallness parameter q, and the harmonic coefficients j2 n - The 
computed values of the six radii for each geoid (rotation period) follow from equation (3) , and the offset 
along z is then added to the calculated radii. 

The residuals in radii at 100 mbar (the measured occultation radii minus the computed radii derived from 
the geoid) are collected into a six dimensional column matrix y, and the linear system y = Ax is inverted to 
obtain corrections to the assumed values of the three parameters. The inverse of the 6x3 rectangular matrix 
A is the pseudo inverse in the singular value decomposition (SVD) method (Lawson & Hanson, 1974). The 
process is repeated to convergence. The matrix A is the design matrix in the least squares process. It is 
computed by applying small finite differences to the three parameters and by evaluating the linear effect on 
the residuals. The resulting A matrix is well conditioned and the least squares solution is full rank three. 
The converged covariance matrix is (A T A)~ 1 af, with T indicating a transpose. The variance of is estimated 
by (y T y)/3 with the result ay = 2.38 km. Using this approach one does not need to set the polar radius a 
priori. Once the solution for the equatorial radius is found, the value of the polar radius b is obtained from 
equation (3) with cos9 set to one. 

A best fit to the occultation radii is found for a rotation period of 9 h 55 m 27 s ± 15 s, a rotation period 
in excellent agreement with Jupiter's measured rotation period. The equatorial radius a, polar radius b, and 
the offset AZ for this rotation period are found to be 71,540.8 km ± 1.4 km, 66,895.2 km ± 3.9 km, and 
5.9 km ±1.7 km, respectively. It should be noted that the errors given are formal mathematical estimates 
and that they do not reflect the real uncertainty in the rotation period estimate which can be of the order 
of about 1 minute. 

The planetocentric latitude 4> an d the radii r of the 100 mbar isosurface are listed in Table 2. We also 
do an extrapolation of the radii to the one-bar level, for purposes of comparing with the elevations inferred 
from the wind data to the 1 bar level. This extrapolation is accomplished by subtracting 47.82 km along 
the vertical to the geoid and projecting it on the geocentric radius. The estimate of the 47.82 km height of 
the 100 mbar level above the one-bar level in the atmosphere is obtained by interpolating in Table 8.2 in 
Lodders & Fegley (1998). We get fairly good agreement with the 1 bar isosurface (1 bar shape) derived in 
the previous section based on the wind data. The fourth column ('Residuals') presents the height difference 
between the calculated geoid for the 100 mbar pressure level (with the offset AZ included) and the measured 
occultation radii for the best fit model (P = 9h 55m 27s) . For our best fit rotation period, the residuals do 
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not exceed 3 km. As in previous sections, but using another source for Jupiter's physical shape, we find that 
the rotation period of Jupiter's deep interior can be ascertained by minimizing the dynamical heights of the 
(100 mbar) isosurface. 

Both radio occultation radii and wind radii for Saturn have been fit in paper I for four rotation periods 
from lOh 32m 35s to lOh 41m 35s. The approach used here for Jupiter can be applied to the Saturn 
occultation data. As for Jupiter, this calculation uses a physical shape which is independent of wind data. 
Unfortunately, the Pioneer and Voyager occultation radii arc not well distributed, with four of the five 
latitudes being in the southern hemisphere. As far as we know, there are no published radii from Cassini 
radio occultations. For Saturn we set the Z offset at -5.75 km and search for the rotation period and equatorial 
radius that can best reconstruct the occultation radii. The rotation which best fits the occultation radii is 
found to be lOh 32m 35s ± 17s, essentially, the same period as found in paper I. The equatorial radius a 
and polar radius b for this rotation period are found to be 60,355.4 km ± 2.9 km and 54,438 km ± 10 km, 
respectively. The five occultation radii and their residuals for the best fit model (P = lOh 32m 35s) are 
shown in Table 3. The larger residual values for Saturn, as compared to Jupiter, may be a result of more 
dynamic horizontal wind motions for Saturn, and hence a larger deviation of the isobaric surface from the 
geoid surface. As a result, the uncertainty in Saturn's rotation period might be larger than that of Jupiter 
and could be about a couple of minutes. 



3.1 Jupiter and Saturn Rotation Periods from Equilibrium Theory 



If one uses the measured (independent of rotation period) polar and equatorial radii of a planet, the rotation 
period which best reproduces the planetary shape can be found using equilibrium theory (Zharkov & Tru- 
bitsyn, 1978). The measured radii of Jupiter and Saturn are given in table 4. This treatment of the effective 
potential assumes a negligible departure from solid-body rotation within the atmosphere and interior and 
therefore assumes that the effect of the winds on the planetary shape is minimal. In principle, accounting 
for differential rotation can influence the planetary figure as well as the gravitational harmonics (Hubbard, 
1982; 1999). However, as we will show below, the shapes of Jupiter and Saturn are found to be accurately 
reproduced using a solid-body rotation. 

The effective potential of the planet (equation (2)) can be used to determine the potential at both the 
pole and the equator. By equating the two results for the potential, the rotation rate of the planet can be 
found. To third order the rotation period is given by (Zharkov & Trubitsyn, 1978), 



T rot = 27T 



2GM 




1- Ji 



J 4 



J 6 



5 

16 



-1/2 



(11) 

where r e and r p are the equatorial and polar radii, respectively. 

The equatorial radius used should be consistent with the reference equatorial radius used for the calcu- 
lation of the gravitational coefficients. For Jupiter we use Jacobson's (2003) reference equatorial radius and 
gravitational harmonics as listed in Table 4. Both the equatorial and polar radii correspond to the 1 bar 
pressure level. Jupiter's rotation period is found to be 9h 55m 20s, only 10 seconds shorter than Jupiter's 
measured radio period. This agreement can be taken as confirmation that the dynamical deformation of 
Jupiter's figure by its differential wind motion is negligible compared with the equatorial bulging due to 
the planet's basic rotation. We repeat this exercise using values of the 100 mbar pressure level, for which 
Jupiter's equatorial and polar radii are 71,541 km and 66,896 km, respectively (Lindal et al., 1981). The 
adjusted gravitational coefficients for Jupiter at the 100 mbar pressure level and its radii are listed in table 
4. Using the values for the 100 mbar pressure level a rotation period of 9h 55m 20s is derived, identical to 
the one derived for the 1 bar pressure level. 
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For Saturn we use the gravitational field measured by Cassini ( Jacobson ct al. , 2006) . These gravity data 
however, correspond to a reference equatorial radius of 60,330 km and not to Saturn's equatorial radius of 
60,268 km at the 1 bar pressure level . We therefore adjust the gravitational harmonics to fit an equatorial 
radius of 60,268km and a polar radius of 54,364 km (see table 4). We find that Saturn's shape is reproduced 
for a rotation period of lOh 31m 49s. This rotation period is less than one minute shorter than the rotation 
period derived by the dynamical heights minimization method. For the 100 mbar pressure level we take 
the equatorial and polar radii to be 60,367 km and 54,438 km, respectively (Lindal et al., 1985). The 
corresponding gravitational coefficients are given in table 4. The rotation period which best reproduces 
Saturn's shape at the 100 mbar pressure level is found to be lOh 31m 50s, a period in excellent agreement 
with the value derived for the shape at the 1 bar pressure level. 

4 Summary and Conclusions 

Our work suggests that the rotation period of Jupiter's deep interior can be obtained by minimizing the 
dynamical heights of its wind level pressure isosurface. We find that a rotation period of 9h 54m 29.7s ±lm 
minimizes the dynamical heights of the 1 bar isosurface. A simpler version of this method minimizes the 
dynamical height difference between the equator and the pole (section 2.1). This procedure provides the 
same rotation period. The derived rotation period, however, is based on the dynamical height calculation 
which is valid for a vertically hydrostatic, pressure-balanced flow, and assumes a vertically uniform velocity. 
Measurement errors may lead to discrepancies as well. As a result, the calculation of the dynamical heights 
can be in error by a few kilometers, resulting in an error of about 1 minute in rotation period. Thus, the 
results suggest that Jupiter's rotation period can be derived by minimizing the dynamical heights of its 
isobaric surface, as suggested in paper I for Saturn. 

We investigate another way of examining the dynamical heights minimization by fitting Pioneer and 
Voyager radio occultation radii at the 100 mbar level in the atmosphere. Again, we search for a rotation 
period which minimizes Jupiter's (100 mbar) isosurface, and obtain a best fit rotation period of 9h 55m 27s, 
consistent with Jupiter's IAU system III rotation of 9h 55m 29.7s. We also use Pioneer and Voyager radio 
occultation radii at the 100 mbar level for Saturn, and search for the rotation period that minimizes the 
dynamical heights. A rotation period of lOh 32m 35s is derived, identical to the rotation period derived in 
paper I. 

The successful application of the dynamical height minimization method to Jupiter lends support to its 
relevance for Saturn. We therefore conclude that the rotation period of Saturn's deep interior is about lOh 
32m. Although we conclude that Jupiter and Saturn rotation periods can be derived by minimizing the 
dynamical heights of their pressure isosurfaces, it is not clear whether the same method can be applied to 
Uranus and Neptune. An analysis of the dynamical height minimization method for Uranus and Neptune 
is more challenging due to larger uncertainties in the planetary physical parameters (radii, gravitational 
moments, rotation periods, etc.) as well as lack of accurate wind measurements and published occultation 
radii. 

Yet another approach for deriving the rotation periods of Jupiter and Saturn uses equilibrium theory. 
Rotation periods of 9h 55m 20s and lOh 31m 49s are found for Jupiter and Saturn, respectively. 

The faster rotation of Saturn advocated here is also supported by an independent analysis of the dynamical 
state of Saturn's atmosphere by Read et al. (2009). These authors analyze the latitudinal distribution of 
potential vorticity on Saturn and show that at many latitudes the atmosphere is close to a state of marginal 
stability with respect to Arnol'd's second stability theorem, an atmospheric configuration allowing long 
planetary Rossby waves that reach deep onto the interior to be coherent across adjacent alternating jet 
streams. Read et al. (2009) infer a rotation period of Saturn equal to lOh 34m 13s ± 20s, about one and a 
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half minutes longer than the ~ lOh 32m 35s period proposed here. The Read et al. (2009) analysis provides 
a dynamical basis for the faster rotation of Saturn. 

Finally, we find that the dynamical deformation of Jupiter by differential wind motions is negligible. 
By assuming the same for Saturn, we infer a rotation period of lOh 32m. We suggest that the Jovian and 
Saturnian figures are dominated by the planetary rotation and that the zonal winds have a minor effect on 
the shapes of both planets. 

The successful inference of Jupiter's known rotation period and the agreement in the values of Saturn's 
rotation period, using different approaches, support the conclusion that Saturn's rotation rate is faster than 
previously thought. 
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Jupiter Rotation Period 


9h 57m 29.7s 


9h 55m 29.7s 


9h 54m 29.7s 


9h 53m 29.7s 


a (km) 


71468.7 
0.0885127 


71491.5 
0.0891935 


71503. 
0.089537 


71514.5 
0.0898826 


J 2 (io- 6 ) 


14706.0 


14696.6 


14691.9 


14687.2 


J 4 (io- 6 ) 


-587.906 


-587.157 


-586.8 


-586.4 


J 6 (io- 6 ) 


34.317 


34.3 


34.2 


34.2 



Table 1: Physical parameters for reference geoids with rotation periods between 9h 57m 29.7s and 9h 53m 
29.7s. 



Occultation Radii for Jupiter 



4> (deg) 


r, 100 mbar 


r, one bar 


Residuals 


Spacecraft 




(km) 


(km) 


(km) 




60.3 


67933.93 


67886.18 


-0.94 


Pio 10 


28.0 


70415.08 


70367.34 


2.97 


Pio 10 


19.8 


70943.95 


70896.17 


-0.98 


Pio 11 


0.07 


71538.61 


71490.79 


-2.23 


Voy 1 


-10.1 


71378.73 


71330.92 


1.13 


Voy 1 


-71.8 


67293.64 


67245.85 


0.04 


Voy 2 



Table 2: Jupiter's planetoccntric latitude (f> and radius r obtained by a digitization of Fig. 7 in Lindal et 
al. (1981) as described in the text. The radii at one bar are extrapolated from the measured radii at 100 
mbar under the assumption of no vertical wind gradients. The residuals are for the 100 mbar radii and are 
from a converged fit with the period as a free parameter. The best fit model is found for a rotation period 
of 9h 55m 27s. 



Occultation Radii for Saturn 



4> (deg) 


r, 100 mbar 


r, one bar 


Residuals 


Spacecraft 




(km) 


(km) 


(km) 




30.5 


58545.4 


58452.8 


-4.6 


Voy 2 


-2.4 


60353.5 


60260.9 


11.1 


Voyl 


-9.8 


60138.0 


60045.4 


-2.5 


Pio 11 


-26.6 


58913.4 


58820.8 


-20.4 


Voy 2 


-71.2 


54948.4 


54855.8 


2.7 


Voy 1 



Table 3: Saturn's latitude <j) and radius r. The radii at one bar are extrapolated from the measured radii 
at 100 mbar under the assumption of no vertical wind gradients. The residuals are for the 100 mbar radii 
and are from a converged fit with the period as a free parameter. The best fit model is found for a rotation 
period of lOh 32m 35s. 
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Parameter 



Jupiter 1 bar Jupiter 100 mbar Saturn 1 bar Saturn 100 mbar 



r e (km) 
T p (km) 
GM (km 3 s~ 2 ) 

J 2 (i(r 6 ) 
J 4 (io- 6 ) 
J 6 (io- 6 ) 



71,492 
66,854 
126,712,765 
14696.43 
-587.14 
34.25 



71,541 
66,896 
126,712,765 
14676.31 
-585.53 
34.11 



60,268 
54,364 
37,931,207.7 
16324.2 
-939.69 
86.67 



60,367 
54,438 
37,931,207.7 
16270.75 
-933.54 
85.82 



Table 4: Data are taken from JPL database: 
(2006), and Lindal et al., (1981; 1985). 



http://ssd.jpl. nasa.gov[ Jacobson (2003), Jacobson et al. 
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Figure 1: Saturn's altitudes of the 100 mbar isobaric surface above a reference geoid with three different 
periods (10 hours, 32 min, 35 s (PI); 10 hours, 35min, 35 s (PII); and 10 hours, 38 min, 35 s (PHI)), and 
a polar radius of 54,438 km. The circles represent radii obtained by radio occultation measurements with 
the Pioneer 11, Voyager 1, and Voyager 2 spacecraft (Lindal et al., 1985). The solid curves represent a 100 
mbar isosurface in geostrophic balance based on zonal-wind data obtained from Voyager. The vertical line 
represents a reference geoid unperturbed by zonal winds. 
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Figure 2: Average zonal wind velocity profile for Jupiter, from HST images between 1995 and 1998, from 
Garci'a-Melendo & Sanchez-Lavega (2001). Positive wind velocities are eastward. 
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Figure 3: Jupiter's altitudes of the 1 bar isosurface above a reference geoid for four different rotation periods 
(9 hours, 53 min, 29.7 s (PI), 9 hours, 54 min, 29.7 s (P2), 9 hours, 55 min, 29.7 s (P3), and 9 hours, 57 min, 
29.7 s (P4)), and a polar radius of 66,854 km. The vertical line represents a reference geoid unperturbed by 
zonal winds. 
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